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We have proposed the non-Abehan dual superconductivity picture for quark confinement in the SU(3) Yang-Mills 
(YM) theory, and have given numerical evidences for the restricted-field dominance and the non-Abelian magnetic 
monopole dominance in the string tension by applying a new formulation of the YM theory on a lattice. To establish 
the non-Abelian dual superconductivity picture for quark confinement, we have observed the non-Abelian dual 
Meissner effect in the SU(3) Yang-Mills theory by measuring the chromoelectric flux created by the quark-antiquark 
source, and the non-Abelian magnetic monopole currents induced around the flux. We conclude that the dual 
superconductivity of the SU(3) Yang-Mills theory is strictly the type I and that this type of dual superconductivity is 
reproduced by the restricted field and the non-Abelian magnetic monopole part, in sharp contrast to the SU(2) case: 
the border of type I and type II. 



Xth Quark Confinement and the Hadron Spectrum, 
October 8-12, 2012 

TUM Campus Garching, Munich, Germany 



* Speaker. 



non-Abelian dual superconductivity for SU(3) Yang-Mills theory 



Akihiro Shibata 



1. Introduction 

Quark confinement follows from the area law of the Wilson loop average. The dual superconductivity is the 
promising mechanism for quark confinement In many preceding works, the Abelian projection [||] was used 
to peifoiTn numerical analyses, which exhibited the remarkable results such as Abelian dominance [^, magnetic 
monopole dominance [Q], and center vortex dominance [||] in the string tension. However, these results are obtained 
only in special gauges: the maximal Abelian (MA) gauge and the Laplacian Abelian gauge within the Abelian pro- 
jection, which breaks the gauge symmetry as well as color symmetry (global symmetry). 

In order to overcome the shortcomings of the Abelian projection, we have presented a new lattice formulation 
of SU (N) Yang-Mills (YM) theory in the previous papers ^ (as a lattice version of the continuum formulations 



for SU (2) and [IC] for SU (N)), which gives a decomposition of the gauge link variable suited for extracting the 



dominant modes for quark confinement in the gauge independent way. In the case of SU (2), the decomposition of the 



gauge link variable was given on a lattice [ 1 1 ] as a lattice version of the Cho-Duan-Ge-Faddeev-Niemi decomposition 
For the gauge group G = SU (N) {N > 3), it was found that the extension of the decomposition from SU (2) to 
SU (N) (N > 3) is not unique and that there are a number of possible ways of decompositions discriminated by the 



stability subgroup H of G, while there is the unique option of H = U{1) in the SU (2) case []12[]. 

For the case of G = SU (3), in particular, there are two possibilities which we call the maximal option and the min- 
imal option. The maximal option is obtained for the stability group H = U{1) xU{l), which enables us to give a gauge 



invariant version of the MA gauge as the Abelian projection [13, 15]. The minimal one is obtained for the stability 



group H = U{2) = SU (2) xU{l), which is suited for representing the Wilson loop in the fundamental representation 



as derived from the non-Abelian Stokes theorem [14]. In the static potential for a pair of quark and antiquark in the 



fundamental representation, we have demonstrated in [16] and [|T7|]: (i) the restricted-field dominance or "Abelian" 
dominance (which is a gauge-independent (invariant) extension of the conventionally called Abelian dominance): the 
string tension ay obtained from the decomposed V-field (i.e., restricted field) reproduced the string tension afuu of 
the original YM field, Gy/OMi = 93 ± 16%, (ii) the gauge-independent non-Abelian magnetic monopole dominance: 
the string tension Oy extracted from the restricted field was reproduced by only the (non-Abelian) magnetic monopole 
part amon, Omon/Ov = 94 ±9%. 

In this paper, we give further evidences for establishing the non-Abelian dual superconductivity picture for quark 
confinement in SU (3) Yang-Mills theory claimed in [ [TtI ] by applying our new formulation to the SU (3) YM theory 
on a lattice. First, we study the dual Meissner effect by measuring the distribution of chromo-flux created by a pair of 
static quark and antiquai^k. We compare the chromo-flux of the original Yang-Mills field with that of the restricted field 
and examine if the restricted field corresponding to the stability group H = U{2) reproduces the dual Meissner effect, 
namely, the dominant part of the chromoelectric field strength of SU (3) Yang-Mills theory. Second, we measure the 
possible magnetic monopole current induced around the flux connecting a pair of static quark and antiquark. Third, we 
focus on the type of dual superconductivity, i.e., type I or type II. In the SU (2) case, the extracted field corresponding 
to the stability group H = U{\) reproduces the dual Meissner effect, which gives a gauge invariant version of MA 
gauge in the Abelian projection, as will be given in [jTsI]. In this paper, we find that the dual superconductivity of the 



SU (3) Yang-Mills theory is indeed the type I, in sharp contrast to the SU (2) case: the border of type I and type II [18]. 



2. Lattice formulation 

We focus our studies on confinement of quarks in a specific representation, i.e., the fundamental representation. 
For this purpose, we consider the Wilson loop operator for obtaining the quark potential, magnetic monopole current 
and chromo-field strength in a gauge invariant way. The Wilson loop operator is uniquely defined by giving a repre- 
sentation, to which the source quark belongs. A remarkable fact is that the Wilson loop operator in the fundamental 
representation leads us to the minimal option in the sense that it is exactly rewritten in terms of some of the variables 
(i.e., the color field n and the V field) which are the same as those adopted in the minimal option, as shown in the 
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process of deriving a non-Abelian Stokes theorem for the Wilson loop operator by Kondo [[l4|]. Therefore, we use the 
reformulation of the Yang-Mills theory in the minimal option to calculate the average of the Wilson loop operator in 
the fundamental representation. We give a brief summary of a new formulation of the lattice SU (3) YM theory ^. 

For the original SU(3) gauge link variable Ux,fi G SU (3), we wish to decompose it into new variables Vx,^ and 
Xc,^ which have values in the SU{3) group, i.e., X^^^ G SU{3), Vx,ii £ SU{3), Ux,^ = Xx^^iVx^n G SU{3), so that Vx.^ 
could be the dominant mode for quark confinement, while Xx,^ is the remainder. In this decomposition, we require 
that Vx,ii is transformed in the same way as the original gauge link variable Ux^^i and as a site variable by the full 
SU{3) gauge ti^ansformation Hx'. Ux^^ — > U[ y = D-xUx^^dI^^, 

^x,ii ^ ^v,v — ^xyx,^l^l+jJ,l ^x,ll ^ ^x,v — ^xXx,ll^l- (2.1) 

We introduce the key variable hx called the color field. In the minimal option, the color field is defined by = 
(§(AV2)i§^ G Lie[SU{3)/U{2)], with being the Gell-Mann matrix and ^ the SU{3) group element. Then, the 
decomposition is uniquely determined from Eqs.(p!^, if the color field hv is specified |^: 

X,,^ = det(L,,^)l/^ Vx,^ =Xl^Ux,^ = gxLx,iiUx,^. (2.2a) 

1x,ii = {Lx,^Ll^y^^^ Lx,^, Lx,^ = ^1 + \/f ('^-^ + Ux,^hx+iiUl^i) + ^\ixUx,^^x+iiUl^ ■ (2.2b) 

In order to determine the configuration {h^} of color fields, we use the reduction condition [^, ^ which guarantees 
that the new theory written in terms of new variables {Xx^^,Vx,ii) is equipollent to the original YM theory. Here, we 
use the reduction condition: for a given configuration of the original link variables Ux,fi, color fields {h^} ai^e obtained 
by minimizing the functional: 

^red[{h.}] = £tr{(D^[[/,,^]h,-)^(D^[t/.,^,]h,)} • (2.3) 

x,ll 



3. Method and results 

We generate configurations of the YM gauge link variable {Ux^^ } using the standard Wilson action on a 24^ lattice 
at j8 = 6.2. The gauge link decomposition is obtained according to the framework given in the previous section: the 



color field configuration {hx} is obtained by solving the reduction condition of minimizing the functional eq.(2.3) for 



each gauge configuration {Ux^p}, and then the decomposed variables {Vv,^}> {Xr,^} ai^e obtained by using the formula 



eq.(|2.2|). In the measurement of the Wilson loop average, we apply the APE smearing technique to reduce noises []20|]. 



3.1 Dual Meissner effect 

We investigate the non-Abelian dual Meissner effect as the mechanism of quark confinement. In order to extract 
the chromo-field, we use a gauge-invariant correlation function proposed in [p^]: The chromo-field created by a quark- 
antiquark pair in SU {N) Yang-Mills theory is measured by using a gauge-invaiiant connected con^elator between a 
plaquette and the Wilson loop (see left panel of Fig.|I|): 

f.v[f|:=s-,/^Pw, p.:J"(f^f'f»-l"^f'"r». (3.1) 

where Fpy\U] is the gauge-invariant chromo-field strength, j3 := 2N /g^ the lattice gauge coupling constant, W the 
Wilson loop in Z-T plane representing a pair of quark and antiquark. Up a plaquette variable as the probe operator 
to measure the chromo-field strength at the point P, and L the Wilson line connecting the source W and the probe 
Up. Here L is necessary to guarantee the gauge invariance of the correlator piy and hence the probe is identified with 
LUpL^. The symbol denotes the average of the operator in the space and the ensemble of the configurations. 
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color f[ux: Original Yang-Mills (L/e =8 , z/e =4 ) 



color flux: restricted field (L/e =8 , z/e =4 ) 
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Figure 1: Measurement of components of the chromoelectric field E and chromomagnetic field B as functions of the distance 
y from the z axis, (left) The gauge invariant connected correlator (UpLWL^) between plaquette U and Wilson loop W. (Center 
panel) the original SU (3) YM field, (Right panel) the restricted field. 




Figure 2: The distribution in Y-Z plane of the chromoelectric field E. connecting a pair of quark and antiquark: (Left panel) 
chromoelectric field produced from the original YM field, (Right panel) chromoelectric field produced from the restricted field. 

We measure con^elators between the plaquette Up and the chromo-field strength of the restricted field Vx,^ as well 
as the original YM field JT,,^. See the center panel of Fig. [T[ Here the quark and antiquark source is introduced as 
8x8 Wilson loop (W) in the Z-T plane, and the probe (Up) is set at the center of the Wilson loop and moved along the 
y-direction. The center and right panel of Fig. |] show respectively the results of measurements for the chromoelectric 
and chromomagnetic fields F^y[U] for the original SU (3) field U and ^fiv[V^] for the restricted field V, where the field 
strength f^vi^] is obtained by using Vx,^ in eq( |3.1[ ) instead of Ux,^. We have checked that even if W[?7] is replaced by 
together with replacement of the probe LUpL' by the corresponding V version, the change in the magnitude of 
the field strength fuv remains within at most a few percent. 

From Fig.[l]we find that only the component of the chromoelectric field {Ex,Ey,E^) = (^^10,^201^30) connecting 
q and q has non-zero value for both the restricted field V and the original YM field U. The other components are 
zero consistently within the numerical errors. This means that the chromomagnetic field {Bx,By,B^) = (^23,^31,^12) 
connecting q and q does not exist and that the chromoelectric field is parallel to the z axis on which quark and antiquark 
are located. The magnitude E^ quickly decreases in the distance y away from the Wilson loop. 

To see the profile of the nonvanishing component of the chromoelectric field in detail, we explore the distribu- 
tion of chromoelectric field on the 2-dimensional plane. Fig. ^ shows the distribution of E. component of the chromo- 
electric field, where the quark-antiquark source represented as 9 x 11 Wilson loop W is placed at (F,Z) = (0,0), (0,9), 
and the probe U is displaced on the Y-Z plane at the midpoint of the T-direction. The position of a quark and an an- 
tiquark is marked by the solid (blue) box. The magnitude of E^ is shown by the height of the 3D plot and also the 
contour plot in the bottom plane. The left panel of Fig. ^ shows the plot of for the SU (3) YM field U, and the right 
panel of Fig. ^for the restricted-field V. We find that the magnitude E, is quite uniform for the restricted part V, while 
it is almost uniform for the original part U except for the neighborhoods of the locations of q, q source. This difference 
is due to the contributions from the remaining part X which affects only the short distance, as will be discussed in the 
next section. 



non-Abelian dual superconductivity for SU(3) Yang-Mills theory 



Akihiro Shibata 



magnetic current Jm x 




distance y/e 

Figure 3: The magnetic-monopole current k induced around the flux along the z axis connecting a quark-antiquark pair (Left 
panel) The positional relationship between the chromoelectric field E- and the magnetic current k. (Right panel) The magnitude 
of the chromo-electronic current E, and the magnetic current J,n — |k| as functions of the distance y from the z axis. 



3.2 Magnetic current 

Next, we investigate the relation between the chromoelectric flux and the magnetic current. The magnetic(- 
monopole) current can be calculated as 

k = *dF[\], (3.2) 

where F[\] is the field strength (2-form) of the restricted field (1-form) V, d the exterior derivative and * denotes the 
Hodge dual operation. Note that non-zero magnetic current follows from violation of the Bianchi identity (If the field 
strength was given by the exterior derivative of V field (one-form), F[\] = d\, we would obtain k =*d^\ = 0). 

Fig. m shows the magnetic current measured in X-Y plane at the midpoint of quark and antiquark pair in the Z- 
direction. The left panel of Fig. ^ shows the positional relationship between chromoelectric flux and magnetic current. 
The right panel of Fig. |3| shows the magnitude of the chromoelectric field (left scale) and the magnetic current k 
(right scale). The existence of nonvanishing magnetic current k around the chromoelectric field supports the dual 
picture of the ordinary superconductor exhibiting the electric current J around the magnetic field B. 

In our formulation, it is possible to define a gauge-invariant magnetic-monopole current by using V-field, 
which is obtained from the field strength =^^v[V] of the field V, as suggested from the non-Abelian Stokes theorem 



[]14|]. It should be also noticed that this magnetic-monopole current is a non-Abelian magnetic monopole extracted 
from the V field, which corresponds to the stability group H = U{2). The magnetic-monopole current defined in 
this way can be used to study the magnetic current around the chromoelectric flux tube, instead of the above definition 



of k (3.2). The comparison of two monopole currents k will be done in the forthcoming paper. 



3.3 Type of dual superconductivity 

Moreover, we investigate the QCD vacuum, i.e., type of the dual superconductor. The left panel of Fig.^ is the 
plot for the chromoelectric field E^ as a function of the distance y in units of the lattice spacing e for the original SU (3) 
field and for the restricted field. 

In order to examine the type of the dual superconductivity, we apply the formula for the magnetic field derived 
by Clem [ ^2| ] in the ordinary superconductor based on the Ginzburg-Landau (GL) theory to the chromoelectric field 
in the dual superconductor. In the GL theory, the gauge field A and the scalar field ^ obey simultaneously the GL 
equation and the Ampere equation: 

[d^ - iqA^){d^ - iqA^)(j)+X{r^ " r]') = 0, (3.3) 



Usually, in the dual superconductor of the type II, it is justified to use the asymptotic form Ko{y/X) to fit the 
chromoelectric field in the large y region (as the solution of the Ampere equation in the dual GL theory). However, 
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Figure 4: (Left panel) The plot of the chromoelectric field E- versus the distance y in units of the lattice spacing e and the 



fitting as a function E-Jj) of y according to ( 3.8). The red cross for the original 5t/(3) field and the green square symbol for 



the restricted field. (Right panel) The order parameter reproduced as a function <^{y) of y according to (3.6), togather with the 
chromoelectric field E-Sy). 
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0.598(5) 


1.878(4) 


1.672(14) 


3.14(9) 


3.75(12) 


4.36(30) 
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restricted field 


0.435(3) 


0.547(7) 


1.787(5) 


1.828(23) 


3.26(13) 


3.84(19) 


2.96(30) 


0.48(2) 



Table 1: The properties of the Yang-Mills vacuum as the dual superconductor obtained by fitting the data of chromoelectric field 
with the prediction of the dual Ginzburg-Landau theory. 



it is clear that this solution cannot be applied to the small y region, as is easily seen from the fact that K^iylX) — oo 
as J — )• 0. In order to see the difference between type I and type II, it is crucial to see the relatively small y region. 
Therefore, such a simple form cannot be used to detect the type I dual superconductor. However, this important aspect 



was ignored in the preceding studies except for a work [ |27| | 



On the other hand, Clem [ ]22| ] does not obtain the analytical solution of the GL equation explicitly and use an 
approximated form for the scalar field ^ (given below in ( |3.6| )). This fonii is used to solve the Ampere equation 
exactly to obtain the analytical form for the gauge field and the resulting magnetic field B. This method does not 
change the behavior of the gauge field in the long distance, but it gives a finite value for the gauge field even at the 
origin. Therefore, we can obtain the formula which is vahd for any distance (core radius) y from the axis connecting 
q and q: the profile of chromoelectric field in the dual superconductor is obtained: 

O 1 Kq{R/X) 



^'W = i^i7§if.«=^/?K^. (3.5) 

provided that the scalar field is given by (See the right panel of Fig.Q) 

O 1 y 

2n V2A ^/y^ + 

where Ky is the modified Bessel function of the v-th order, A the parameter con^esponding to the London penetration 
length, a variational parameter for the core radius, and external electric flux. In the dual superconductor, we 
define the GL parameter K as the ratio of the London penetration length A and the coherence length <^ which measures 
the coherence of the magnetic monopole condensate (the dual version of the Cooper pair condensate): K = X/^. It is 
given by ||2|] 

K = V2^^l-Kl,ii;/l)/K^iUX). (3.7) 



According to the formula Eq.(3.5), we estimate the GL parameter K for the dual superconductor of SU (3) YM 
theoty, although this formula is obtained for the ordinary superconductor of U (1) gauge field. By using the fitting 
function: 

Oil ,1 



Eiy)=aMVWfT?). „ = -_^^, (, = -, (3.8) 
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we obtain the result shown in Table |I} The superconductor is type I if k: < Kc, while type II if K" > jq, where the 
critical value of GL parameter dividing the type of the superconductor is given hy Kc = I/V2 ~ 0.707. Our data 
clearly shows that the dual superconductor of SU (3) Yang-Mills theory is type I with 

K- = 0.45 ± 0.01 , A =0.1207 ±0.017 fm, ,^ = 0.2707 ± 0.086 fm, (3.9) 

where we have used the string tension aphys = (440MeV)^, and data of lattice spacing is taken from the TABLE I in 



Ref.[21|. This result is consistent with a quite recent result obtained independently by Cea, Cosmai and Papa [27]. 
Moreover, our result shows that the restricted part plays the dominant role in detemiining the type of the non-Abelian 
dual superconductivity of the SU (3) Yang-Mills theory, i.e., type I with 

k: = 0.48 ± 0.02 , A = 0.132 ±0.03 fm , <^ = 0.277 ±0.014 fm. (3.10) 

This is a novel feature overlooked in the preceding studies. Thus the restricted-field dominance can be seen also 
in the determination of the type of dual superconductivity where the discrepancy is just the normalization of the 
chromoelectric field at the core y = 0, coming from the difference of the total flux <I>. These are gauge-invariant 
results. Note again that this restricted-field and the non-AbeUan magnetic monopole extracted from it reproduce the 



string tension in the static quark-antiquark potential []17[]. 



Our result should be compared with the result obtained by using the Abelian projection: Y. Matsubai^a et. al [25| 



suggests K = 0.5 ~ l(which is j3 dependent), border of type I and type II for both SU (2) and SU (3). In SU (2) case, on 



the other hand, there are other works [ 24 1 which conclude that the type of vacuum is at the border of type I and type 



II. We should mention the work [|^] which concludes that the dual superconductivity of SU (3) Yang -Mills theory is 



type II with K" = 1.2 ~ 1.3. This conclusion seems to contradict our result for SU (3). If the above formula (3.5) is 



applied to the data of [26], we have the same conclusion, namely, the type I with K = 0.47 ~ 0.50. Therefore, the data 



obtained in [26] are consistent with ours. The difference between type I and type II is attributed to the way of fitting 



the data with the formula for the chromo-field. 



4. Summary and outlook 

We have given further numerical evidences for confirming the non-Abelian dual superconductivity for SU (3) YM 
theory proposed in [17]. For this purpose, we have used our new formulation of SU (3) YM theory on a lattice [g 0] 



to extract the restricted field from the original SU (3) YM field, which has played a dominant role in confinement of 
quarks in the fundamental representation, i.e., the restricted-field dominance and the non-Abelian magnetic monopole 



dominance in the string tension, as shown in the previous studies [17]. 

We have focused on the dual Meissner effect and have measured the chromoelectric field connecting a quark and 
an antiquark for both the original YM field and the restricted field. We have observed the dual Meissner effect in SU (3) 
YM theory, i.e., only the chromoelectric field exists and the magnetic-monopole current is induced around the flux 
connecting a quark and an antiquark. Moreover, we have determined the type of non-Abelian dual superconductivity, 
i.e., type I for SU (3) YM theory, which should be compared with the border of type I and II for the SU(2) YM 
theory. These features are reproduced only from the restricted part. These results confirm the non-Abelian dual 
superconductivity picture for quark confinement. 
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